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Hadron deformation from Lattice QCD * 
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We address the issue of hadron deformation within the framework of lattice QCD. For hadrons with spin greater 
than 1/2 the deformation can be determined by evaluating the charge and matter distributions. Deviation of the 
nucleon shape from spherical symmetry is determined by evaluating the quadrupole strength in the transition 
jN — > A(1232), both in the quenched and in the unquenched theory. 

1. Introduction 



Hadronic matrix elements of two- and three- 
current correlators encode detailed information 
on hadron structure such as quark spatial distri- 
butions, charge radii and quadrupole moments. 
In the non-relativistic limit they reduce to the 
square of the wave function and therefore the 
shape of the hadron can be extracted. The is- 
sue of a deformation in the nucleon and A was 
first examined in the context of the constituent 
quark model £Q with the one gluon exchange giv- 
ing rise to a colour magnetic dipole-dipole inter- 
action. The tensor component of this interac- 
tion causes a D-wave admixture in the nucleon 
and A wave function and leads to a deforma- 
tion. Other mechanisms have been invoked to 
explain hadron deformation: In cloudy bag mod- 
els the deformation is caused by the asymmetric 
pion cloud In soliton models it is thought 
to be due to the non-linear pion field interaction. 
In the constituent quark model it was recently 
proposed to be due to two-body contributions to 
the electromagnetic current arising from the elim- 
ination of gluonic and quark-antiquark degrees of 
freedom Because of the strong experimen- 
tal indications that the nucleon is deformed, it is 
interesting to address the general issue of hadron 
deformation within lattice QCD and try to under- 
stand its physical origin. In the first part of this 
talk we present results on the charge and mat- 
ter density distributions evaluated on the lattice, 
both in the quenched and in the unquenched the- 
ory. A comparison between the charge and mat- 



ter distributions suggests that deformation can 
arise from the relativistic motion of the quarks. 

For the nucleon, being a spin 1/2 particle, 
the spectroscopic quadrupole moment averages to 
zero and therefore it cannot be studied via the 
density distributions. For this reason we look for 
a quadrupole strength in the ~/N — > A transition 
as is done in experimental studies. Spin-parity se- 
lection rules allow a magnetic dipole, Ml, an elec- 
tric quadrupole, E2, or a Coulomb quadrupole, 
C2, amplitude for this transition. If both the nu- 
cleon and the A are spherical, then E2 and C2 
should be zero. Although Ml is indeed the dom- 
inant amplitude, there is mounting experimen- 
tal evidence over a range of momentum transfers 
that E2 and C2 are non-zero 0]. State-of-the- 
art lattice QCD calculations can yield model in- 
dependent results on these form factors and pro- 
vide direct comparison with experiment. The sec- 
ond part of this talk discusses the evaluation of 
the Sachs transition form factors Gmi, Ge2 and 
Qc2 describing -fN — > A. 

2. Hadron wave functions 

For a meson we evaluate the two-current cor- 
relator shown schematically in Fig. ^ Taking the 
current insertions at equal times t\ =t%=t this 
correlator is given by 

C r (r,t) = J d 3 r' (h\ff(r' + r,t)j*(r>,t)\h) (1) 

for a hadron state \h). The current operator is 
given by the normal ordered product 



*This research is supported in part by funds provided by 
the Levendis Foundation 



jr(r, t) =: u(r, t)Tu(r, t) : where T = 



7o 
1 



2 



for the charge and matter correlators respectively. 
We note that no gauge ambiguity arises in this 
determination of hadron wave functions, unlike 
Bethe-Salpeter amplitudes. In the case of baryons 
two relative distances are involved and three cur- 
rent insertions are required as shown in Fig ^ 
However we may consider integrating over one rel- 
ative distance to obtain the one-particle density 
distribution that involves two current insertions 
and is thus evaluated via Eq. ^for t\ = t<x = t. 




, (r'+r,t) 



Meson 



Baryon 



Figure 1. Two- and three- current correlators for 
a meson (left) and a baryon (upper right); Lower 
right: two-current correlator for a baryon after 
integration over one relative distance. These di- 
agrams are drawn assuming anti-periodic bound- 
ary conditions in the time direction. 

All the results on the current-current corre- 
lators have been obtained on lattices of size 
16 3 x 32. We analyse, for the quenched case, 
220 NERSC configurations at j3 = 6.0 and, for 
the unquenched, 150 SESAM configurations |S] 
for k = 0.156 and 200 for k = 0.157 simulated at 
/3 = 5.6 with two dynamical quark flavours. The 
physical volume of the quenched and unquenched 
lattices is approximately the same. 

In Eq. n the current couples to the quark at 
fixed time separation t from the source, which 
must be large enough to sufficiently isolate the 
hadron ground state. We checked, by compar- 
ing data at different t, that excited states are 
sufficiently suppressed at maximal time separa- 
tion t = T/4 in the current-current correlators [7] 
where we have taken anti-periodic boundary con- 
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Figure 2. Quenched charge and matter distribu- 
tions, C 7o (r) and C/(r), at n = 0.153. Left for 
the pion and the rho; Right for the nucleon and 
the A. 

ditions in the temporal direction. This allows 
us to analyse standard full QCD configurations 
that employ anti - periodic boundary conditions. 
Since, for our parameters, local sources produce 
the same results for Cr(r) as smeared ones and 
carry no gauge noise, we use them in this study. 

The quenched charge and matter density dis- 
tributions were studied as a function of the naive 
quark mass, m q = 1/2(1/ k — X/k c ), where k c is 
the hopping parameter value at which the pion 
becomes massless. The dependence of both distri- 
butions on m q is more pronounced for the mesons 
than for the baryons. In particular the pion 
charge and matter distributions broaden as the 
quarks become lighter [817] . In the case of the rho 
the charge distribution shows a stronger depen- 
dence on the quark mass than the matter distri- 
bution, for which no further mass dependence is 
observed for m q ~ 200 MeV. For the nucleon and 
the A + essentially no variation is seen over the 
range of naive quark masses ~ 300 — 100 MeV in- 
vestigated here. The quenched charge and matter 
distributions are compared in Fig.|2at k — 0.153. 
The matter distribution is for all hadrons less 
broad than the charge distribution. The root 
mean square radius (rms) can be extracted by 
evaluating 



E r r 2 Cr(r,t) 
Er^rM) 



(2) 



3 



For non-relativistic states the charge rms radius 
can be written in terms of the form factors as [5] 

(r 2 ) = -^V| 2 ^ n (q)F„i(-q)| q2=0 

n 

= r 2 u + rl -2^dL-d^ (3) 

n 

where the non-relativistic form factors F^(q) = 
(h, q\j^' d \n, 0) are given by 



l"< d q 2 /6 +0(q 4 ) n = h 
q.d^ + 0(q 2 ) n^h 



(4) 



for u- and d- quarks. The dipole-dipole term 
appearing in Eq. [3] can be evaluated if, for in- 
stance, we assume a non-relativistic two-body 
system with equal masses for the u and d quarks 
in the center of mass frame, since then r„ = —Yd 



and < r 2 >= 2(r 



i). Only if we allow current 



insertions at unequal times t\ and t% and take 
<2 — ti large enough so that intermediate excited 
states are sufficiently suppressed then the charge 
radius, for instance of the pion, is obtained from 



3C 7o (q) 



dq 2 



q 2 =0 



dq 2 



q 2 =0 



-g<ri> -(5) 




Figure 3. C r (0,0,z) (Z-curve); C r (x,0,0) and 
C* r (0,y,0) (X,F-curves). 

The rho charge density distribution 8 lOj pro- 
duces a non-zero quadrupole moment or equiva- 
lently a, z—x asymmetry as shown in Fig-EI where 
the z-axis is taken along the spin of the rho. This 
yields a deformation 6 = § (3z ^ r2> - 0.03 ± 0.01 



too 

0.75 
0.50 
! 0.25 



0.00 If 

0.75 
0.50 
0.25 
0.00 



*; — — 1 — 

^ charge 
- ^ R=0.153 - 

\ 

2 v L=0 
5 *• f 
L=2 


\ matter 
- \ K=0.153 - 

\_ 

\ 

* L=0 
\ 

\. 

L=2 


• I /c = 0.155 

^ L = . 


H K = 0.155 
- 1 

\t 

f L=0 
'1 

III H 



0.5 0.0 

r (fm) 



0.5 1.0 

r (fm) 



Figure 4. Decomposition into L = and L = 2 
contributions at k = 0.153 and k = 0.155. Left: 
charge distribution, right: matter distribution. 

in the chiral limit. An angular decomposition of 
the wave function, as shown in Fig. corrobo- 
rates a non-zero charge deformation by producing 
a non-zero L = 2 component ^Ij • No such defor- 
mation is seen for the matter density distribution 
in Figs. El and E|. Since the matter and charge op- 
erators have the same non-relativistic limit, this 
suggests that hadron charge deformation is a rela- 
tivistic effect. This result has strong implications 
for the validity of various phenomenological mod- 
els used in the study of nucleon deformation. 
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Figure 5. Quenched and unquenched charge den- 
sity distributions at K — 0.153 and k = 0.156 
respectively. Errors bars are omitted for clarity. 




Unquenched results obtained at n = 0.156 



(m w /m p = 0.83) and k — 0.157 {m^/r 



0.76) 



can be compared to quenched ones at n = 0.153 
(m^/nip = 0.84) and n = 0.154 (m n /m p = 0.78) 
having similar pion to rho mass ratios. In Figs. [3] 
and El unquenched results for C 7o (r) and C/(r) at 
k = 0.156 are compared to quenched results at 
k = 0.153. The unquenched charge and matter 
distributions show an increase at short distances 
in the case of the pion and the rho whereas for the 
baryons no significant changes are seen. Like in 
the quenched case the unquenched matter distri- 
bution shows a faster fall-off as compared to that 
observed for the charge distribution. Whereas un- 
qucnching leads to an increase in the rho charge 
asymmetry and to a small deformation for the 
A + , as shown by the 3-dimensional contour plots 
of Fig. [7| it has no effect on the matter density dis- 
tribution. This again suggests a relativistic origin 
for the rho charge deformation. Pion cloud con- 
tributions to hadron deformation are expected to 
become significant as we approach the physical 
pion mass and a lattice calculation with lighter pi- 
ons can test if the charge rho asymmetry shows a 
significant increase, as suggested in some models, 
and at the same time if a matter density asym- 
metry shows up. 

3. 7^ — » A transition form factors 

We evaluate the 3-point function, 
{G^" N (t 2 , t\\ p ', p; T)) shown schematically in 
Fig.il as well as (G^'" A (t 2 ,ti;p',p;r)) for the 




Figure 7. Three-dimensional contour plot of the 
charge correlator (red): upper for the rho state 
(cigar shape) and lower for the A + (slightly 
oblate) for two dynamical quarks at K = 0.156. 
Values of the correlator (0.5 for the rho, 0.8 for 
the A + ) were chosen to show large distances but 
avoid finite-size effects. We have included for 
comparison the contour of a sphere (green) . 

reverse A — > jN transition. Exponential decays 
and normalization constants cancel in the ratio 
R a (h,ti;p',p = 

(G^ N (t 2l tv, p P ; r)) (G^ A (t 2) h; p, p '; rt)> 



1 1/2 



(<%GAA(* 2 ,p';r 4 )) (G^(i 2 ,- P; r 4 )> 

n CT (p',p;r;/i), (6) 



t2-tl>l,*l>l 



where G NN and G^ A are the nucleoli and A two- 
point functions evaluated in the standard way. 
The Sachs form factors are obtained by appropri- 
ate combinations of the A spin-index a, current 
direction fi and projection matrices T. For in- 
stance, in the A rest frame p ' = , q = (q, 0, 0) 
HUH) we have 

Qmi = A T^a n 2 (0,-q;+z T 4 ;3) 



Ge2 =A 



1 

3|qf 



n 3 (0,-q;r 1 ;3)+n 1 (0,-q;r 3 ;3) 
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A(p') N(p) 



Figure 8. 7 TV — + A matrix element. The photon 
couples to a quark in the nucleon at time separa- 
tion t\ from the source to produce a A. 

g C2 = ^ : ^n 1 (0,-q;-ir 1 ;4) (7) 

q 

where A is a kincmatical factor. 




Figure 9. (a) Fixed operator and (b) fixed source 
sequential propagators. 



We use two methods to compute the sequen- 
tial propagator needed to build the 3-point func- 
tion: (a) We evaluate the quark line with the pho- 
ton insertion, shown schematically in Fig.|5Ja) by 
computing the sequential propagator at fixed mo- 
mentum transfer q and fixed time t\. We look for 
a plateau by varying the sink-source separation 
time t%. The final and initial states can be chosen 
at the end. (6) We evaluate the backward sequen- 
tial propagator shown schematically in Fig. OJb) 
by fixing the initial and final states. £2 is fixed 
and a plateau is searched for by varying t\ . Since 
the momentum transfer is specified only at the 
end, the 7iV — > A form factors can be evaluated 
at all lattice momenta. 

The parameters of our lattices are given in Ta- 
ble ^ where we used the nucleon mass in the 
chiral limit to convert to physical units, and 
Q 2 = —q 2 is evaluated in the rest frame of the 
A. 



Table 1 



(GeV") 


K 


m n /m p 


Number of confs 


Quenched (5 = 6.0 16 J x 32 


q 




= 0.64 GeV" 


0.64 


0.1530 


0.84 






100 


0.64 


0.1540 


0.78 






100 


0.64 


0.1550 


0.70 






100 


Quenched /3 = 6.0 32^ x 64 


q 




= 0.64 GeV" 


0.64 


0.1550 


0.69 


100 


Quenched j3 = 6.0 32^ x 64 


q 




= 0.16 GeV" 


0.16 


0.1554 


0.64 






100 


0.15 


0.1558 


0.59 






100 


0.13 


0.1562 


0.50 






100 


Unquenched (3 = 5.6 16 J x 32 


6 


q" = 0.54 GeV" 


0.54 


0.1560 


0.83 






196 


0.54 


0.1565 


0.81 






200 


0.54 


0.1570 


0.76 






201 


0.54 


0.1575 


0.68 






200 



We check for finite volume effects by compar- 
ing results in the quenched theory on lattices of 
size 16 3 x 32 and 32 3 x 64 at the same momentum 
transfer at k = 0.1550. Assuming a 1/volume de- 
pendence we find that on the small volumes there 
is a (10 — 15)% correction as compared to the in- 
finite volume result, whereas on the large lattice 
the volume correction is negligible. In Figs. 1101 
and llll we show quenched and unquenched results 
for Qmx an d Qe2 at the same momentum transfer. 
Unquenching tends to decrease Qui and Qe2 but 
leaves the ratio Rem = —QeiIQmi largely un- 
affected for the SESAM quark masses studied in 
this work, giving values in the range of —(2 — 4)%. 
The fact that no increase of Rem is observed 
means that pion contributions to this ratio are 
small for the SESAM pion masses. In Fig. 1121 we 
also show results for Qci with A static for the 
large quenched lattice for which we obtain the 
best signal. Although Qci is within one standard 
deviation of zero, it is positive at all n- values giv- 
ing a negative ratio Rsm = — |q|/2mA QciIQmi 
in agreement with experiment. 

Chiral extrapolation of the results is done lin- 
early in the pion mass squared, since with the 
nucleon or the A carrying a finite momentum, 
chiral logs are expected to be suppressed. The 
values obtained are given in Table |21 and are in 
reasonable agreement with the experimental val- 
ues Rem = -2.1 ± 0.2 ± 2.0 at Q 2 = 0.126 GeV 2 
and Rem = -1.6 ± 0.4 ± 0.4 at Q 2 = 0.52 GeV 2 . 
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quenched 
k=0.153 



k=0.154 



k=0.155 



unquenched 
/c=0.1560 



,a k=0.1570 



till 



k=0.1575 



t,/o 



Figure 10. £ M i at q = (27r/16a, 0, 0) for the 
quenched (left) and for the unquenched (right) 
theory. The photon couples to the quark at 
ti/a = 6 as shown by the arrow. The lines 
show the plateau fit range and bounds obtained 
by jackknife analysis. 

Using the fixed sink method for the sequential 
propagator we obtain the ^-dependence of the 
form factors. In Fig. 1131 we show preliminary re- 
sults for 



c* - 1 



i 



/Ml 



(8) 



1 



(mjv+mA) 2 



obtained using 50 quenched configurations at k = 
0.1554 and 0.1558, and 25 at k = 0.1562. Here 
we have used the ratio 



<%G£ A (; 2 ;p';r 4 )) 



(G NN (2t 1 ;p;Ti)) 
<«5 lJ G^ /i (2t 1 ;p';r 4 )> 



-,1/2 



(9) 



to extract the form factors, since using the sym- 
metric combination given in Eq. would require 
an additional sequential inversion. At the low- 
est momentum transfer we confirm that the re- 
sults from Eqs. El andEI are within error bars, as 
shown by the data at Q 2 = 0.13 GeV 2 included in 
Fig. ED The quark mass dependence is weak es- 
pecially at the higher momenta and a global fit to 
the lattice data may be performed. The fit shows 
that the lattice data are consistent with a simple 
exponential dependence on Q 2 . This fit produces 
a value at Q 2 — in agreement with the Particle 
Data Group result |T3]. 
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— I 1 1— 



k=0.154 



-A A A. 



/c=0.155 



unquenched 

k=0.1560 



.* a A A 

/c=0.1570 



..sill 



k=0.1575 



O 4 8 12 O 4 8 12 16 

V<3 t 2 /a 

Figure 11. Q E 2 at q = (2n/lQa, 0, 0). The nota- 
tion is the same as that of Fig. 1101 



quenched 



I { 



£11 



-^{ : 



, * * i 1 



/c = 0.1562 




'0 2-1 6 8 10 12 11 16 18 

t 2 /a 

Figure 12. Q C 2 at q = (27r/32a, 0, 0) for the 
quenched theory. The photon couples to the 
quark at ti/a — 8 as shown by the arrow. The 
rest of the notation is the same as that of Fie . 1101 

4. Conclusions 

Lattice techniques are shown to be suitable for 
the evaluation of hadron charge and matter den- 
sity distributions. Differences between quenched 
and unquenched results are shown to be insignifi- 
cant, for both the charge and matter density dis- 
tributions. The charge density distribution is, in 
all cases, broader than the matter density. For 
baryons, the lattice indicates a charge radius, 
which is ~ 20% larger than the matter radius. 
The deformation seen in the rho charge distribu- 
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fixed operator (100 confs) 
m n =0.51 GeV (50 confs) 
rn^O.45 GeV (50 confs) 

m =0.37 GeV (25 confs) 



* > I 



1.0 



1.5 



Q 2 (GeV 2 ) 

Figure 13. Q* M1 versus Q 2 in the A rest frame at 
k = 0.1554, 0.1558 and 0.1562. The dashed line 
is an exponential fit to all the lattice data. The 
filled square is the value of Q\ n at Q 2 — 0.0 given 
by the Particle Data Group. 

Table 2 



GeV" 


Qmi (e/2m N ) 


Qe2 (e/2m N ) 


Rem% 


Quenched QCD 


0.64 


1.72(6) 


0.099(19) 


-5.1(1.1) 


0.13 


2.51(6) 


0.104(12) 


-4.5(1.4) 


Unquenched QCD 


0.53 


1.30(4) 


0.050(31) 


-2.8(1.6) 



tion is absent in the matter distribution, both in 
the quenched and the unquenched theory. This 
observation suggests a relativistic origin for the 
deformation and deserves a more extended study, 
with lighter quarks and larger volumes. Nucleon 
deformation is signaled by a non-zero quadrupole 
strength in the transition 7 A — > A. We have cal- 
culated, for the first time using lattice QCD, the 
electric quadrupole form factor to high enough ac- 
curacy to exclude a zero value. The ratio Rem is 
evaluated in the kinematical regime explored by 
experiments. The values extracted in the chiral 
limit are in good agreement with recent measure- 
ments. Although large statistical and systematic 
errors prevent an accurate determination of the 
Coulomb quadrupole form factor and a zero value 
cannot be excluded, our results support a nega- 
tive value of the ratio Rsm in agreement with ex- 
periment. We have shown preliminary results on 
the Q 2 dependence of the magnetic dipole transi- 
tion form factor evaluated with ~ 10% accuracy 
in the regime explored by Jefferson Lab. The 
main conclusion from the analysis of the SESAM 



lattices [S] is that, for pions in the range of 800- 
500 MeV, no unquenching effects can be estab- 
lished for Rem within our statistics. We plan to 
repeat this calculation for lighter pions to inves- 
tigate the pion cloud contributions on the transi- 
tion form factors. 
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